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Abstract. Wc study 1+1 field-generalizations of the rational and elliptic Gaudin models. 
For s\{N) case we introduce equations of motion and L-A pair with spectral parameter on 
^ . the Riemann sphere and elliptic curve. In sl(2) case we study the equations in detail and find 

' the corresponding Hamiltonian densities. The n-site model describes n interacting Landau- 

Lifshitz models of magnets. The interaction depends on position of the sites (marked points 
^ ' on the curve). We also analyze the 2-site case in its own right and describe its relation to 

the principal chiral model. We emphasize that 1+1 version impose a restriction on a choice 
' of flows on the level of the corresponding 0+1 classical mechanics. 
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1 Introduction 

Gaudin model (or Gaudin magnet) was introduced by M. Gaudin [1] as a quasiclassical limit 
I of spin-1/2 chain and was studied via the Bethe ansatz [2l [3]. Let us start with a general 

■ rational model underlying Gaudin magnets. The classical rational Gaudin model is defined by 

. the following quadratic Hamiltonians: 



^a = -^V^, a = l,...,n, (1.1) 



where S"' € sl(2, C), {zi, . . . , Zn} S CP"*^ are marked points and ( ) denotes the trace. 

From the point of view of the Lax pair the model is described by a general Lax matrix which 
is a sl(A^, C)-valued function L{z) on CP"^\{2:i, . . . , Zn} with simple poles at {zi, . . . , Zn} and 
some given residues ReSz^L{z) = 5" E sl(A^, C): 

The generating function of the Hamiltonians is 

a=l aytb 



The first sum in (|1.3p shows that the eigenvalues of S"" are the constant C-numbers. Thus, 
the phase space is a direct product^ of the coadjoint orbits by SL(A^, C) action: M = O^x 



^In fact, there is coadjoint action of SL(A'^, C) on M wliich provides the constraint '^^S'^ — with some 
fixation of SL(A'^, C) action. Then one can make a reduction M — >■ M/ /SL{N, C). But we do not go into details 
of this reduction here. In [4] [5] the examples of the reduction for the Painleve VI equation are discussed. The 
r-matrix of the reduced models satisfies the reflection equations. Thus, the models live on the boundaries of the 
finite lattices. 
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• • • X O". This phase space is naturally equipped with a linear Poisson~Lie structure: 

{5S,5^} = 5'^^^C2^5,^ (1.4) 

7 

where are coefficients in some basis {Tq}: = ^ S^Ta and C^^ are the structure constants 

a 

of sl(A^, C) in this basis. The natural basis is described in the appendix. The Hamiltonians (II. ip 
in sl(A^, C) case are replaced by 

Ha = -^Y.^^^ a = l,...,n. (1.5) 

N ^ Za- Zc 

The dynamics with respect to the Hamiltonians (jl.Sp is given by the following equation^: 

\aa obi 

dt^S" = {H^, 5"} = for a / b. 

These equations of motion can be represented in the Lax form 

dt^L = [L,Ma] (1.6) 
with the Lax pair 

" cc qa 

L{z) = J2^, Ma{z) = . 



In such a generality the model was studied many times. For example, the non-autonomous 
version corresponds to the Schlesinger system of the isomonodromic deformations on a sphere. 
It was studied a hundred years ago [6] . 

In the elliptic case [7j the Lax matrix (jl.2p is replaced by 

n 

a=l a 

where 2; S S,- is a coordinate on an elliptic curve S,- with moduli r. Basis {Tq,} and the 
corresponding Poisson structure is defined in ()B.4p . Functions fai^ ~ ^a) (|B.10p form a basis 
in r(Endy, St-) with a simple pole at Za for some fixed holomorphic vector bundle V of degree 
one. The Poisson structure ()1.4p for the structure constants ()B.7p is related to the existence of 
the r-matrix of the Belavin-Drinfel'd type [8]. The quadratic Poisson structure can be defined 
by the same r-matrix [9]. 

Most of problems natural for integrable systems have been studied for the Gaudin model 
as well. Among them the separation of variables [10], relations to monodromy preserving and 
Knizhnik-Zamolodchikov equations [11], quantum quadratic algebras and bihamiltonian struc- 
tures [12] , time-discrete versions [13] , quantization [14] and Langlands duality [15] . It should be 
mentioned that the elliptic Gaudin model was originally defined by B. Enriquez and V. Rub- 
tsov |16j as an example of the Hitchin-type system [T7]. "Dynamical" case was considered first 
by A. Gorsky and N. Nekrasov [18]. That case corresponded to degree zero vector bundle V 
(that is to nontrivial moduli space of bundles) or to the "spin" extensions of the Calogero model. 



Here we imply some choice of the normalization by the Killing form ( ), see also (|B.6 
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In [19] it was shown that the top-like models and Calogero-type models are related by means of 
the modification procedure (the later changes the degree of V). In this respect, the models are 
equivalent. 

Consideration of particular cases and different types of reductions leads to relations between 
Gaudin model and a number of known integrable systems such as interacting tops [20j , Painle- 
ve VI equation and Zhukovsky-Volterra gyrostat [4], Neumann system |21] . 

With the advent of the inverse scattering method the Lax equations or the zero-curvature 
equations [22] (with spectral parameter) became a main tool for investigation of nonlinear equa- 
tions |231l24j. Different applications and classifications can be found in [25]. In this paper we are 
predominantly interested in the Landau-Lifshitz equation [26] (which describes the continuous 
limit of the XYZ model [271 EH]) and the principal chiral model [29113(1] . 

In [19] a general scheme was suggested for constructing 1 + 1 (or field) generalizations of the 
Gaudin- type models as typical examples of the Hitchin systems. As a by-product of this work 
the field generalization of the elliptic Calogero model was obtainec^ and its equivalence to the 
Landau-Lifshitz equation was shown in terms of the special singular gauge transformations. 

The purpose of the paper is to present explicit L-A pairs for 1 + 1 Gaudin model, to pro- 
pose corresponding Hamiltonian description and to find out relationships between the obtained 
equations and some known models such as the Heisenberg Model, the Landau-Lifshitz equation 
and the principal chiral model. 

In 1+1 models the Lax equations (|1.6p are replaced by the zero-curvature (or Zakharov- 
Shabat) equations: 

dt,L-d,Ma = [L,Ma], 

where L and Ma do not coincide (in general) with those from (II. 6p . It was shown in |19] how to 
construct 1+1 version of L-operator. In particular, L keeps the same form as in (0+1) version 
of the Hitchin systems corresponding to holomorphic vector bundles of degree 1. This class 
of systems is under our consideration in this paper. A general scheme [32] allows to obtain 
densities of the conserved quantities (Hamiltonians) . However, there is a technical problem 
of finding corresponding M-operators. Unfortunately, there is no practical way to get them 
explicitly. For example, in |19j the nontrivial M-operator for the field version of Calogero model 
was obtained by some ansatz. In the same manner M-operators were obtained in [271 [30] for 
the Landau-Lifshitz and the Principal Chiral Models correspondingly. The inverse problem (to 
find mechanical L and M from known field versions) is an easy task - one should put to zero all 
derivatives with respect to the loop variable x. In this respect, there is a correspondence between 
field flows and some choice of flows (M-operators) on the level of classical mechanics. It will be 
shown that the first flows of 1+1 Gaudin hierarchy correspond to "conventional" description of 
flows in the Gaudin mechanics while the second flows arise naturally from some "reformulated" 
version. The later appears as some linear combination of the "conventional" Gaudin flows. 

The paper is organized as follows: in Section [2] we give a standard description of the Gaudin 
model and its flows Ma (Proposition 12. ip . Then the "reformulated" version is suggested in 
the form of linear combinations of {Ma} (Proposition 12. 2^ . In Section [3] we discuss the field 
generalization and find the first (Proposition 13. ip and the second (Proposition 13. 2p flows of 
the 1+1 Gaudin hierarchy. Among other things, we consider a special case of the first flows 
corresponding to the principal chiral model in detail. In Section [J] sl(2, C) case is considered 
(rational - Subsection 14.11 and elliptic - Subsection 14. 2p and the Hamiltonian description is 
obtained. First, we get general formulae for the densities of Hamiltonians via local decomposition 
for the first (Lemma 14. ip and the second (Lemma 14. 2p flows. Secondly, we evaluate these 
densities for 1+1 rational and elliptic sl(2,C) Gaudin model and reproduce previously obtained 
equations of motion (Theorem 14. ip . 

^This result was first obtained by I. Krichever in [31]. 



4 



A.V. Zotov 



The results of the paper can be briefly summarized as follows: 



(0+1) mechanics: 
Gaudin flows {Ha} 



Gaudin flows {Ha} 

( "reformulated version" ) 



(1+1) field version: type of models: 



I''* flows {-Ha,!} 

flows {'Ha,2} 



n-site generalization of 
principal chiral model 

interacting models of 
Landau-Lifshitz type 



The first flows are described by the following equations: 

dt^s'' = [s\0Usi]. 

In "2-site" case and rational limit these are the equations of the principal chiral model: 



dth - kdJo H 



[^i,^o] =0, 



with Iq = S^ + and h = S^-S^. 

The equations for the second flows are of the form (here we put sl(2, C) case and sl(iV, C) is 
considered below): 

di^S'^ - kd^T^- = [S^p(S'^)] + J]K,(^e,(5^)] - [F,a{S'),S% 

[Mvl,S'] + [S\FUSn], (1-8) 



dr S'' 



^[S*", S^] + Yl ^ac{S'^). Note that in "1-site" case n = 1 the first one equation 



where r/" 

in (|1.8p is the Landau-Lifshitz equation (for ti = t): 
k^ 



dtS + ^[S,S, 



[S,p{S)]. 



2 sl(Ar, C) elliptic Gaudin model 
2.1 Standard description 

The phase space of the Gaudin model is a direct product of orbits Oi x • • • x On by the coadjoint 
action of SL(A^, C)). The coordinates {S^} on each orbit S'^ G Oc are chosen to be dual to the 
basis {Tq} of the Lie algebra sl(A^, C). The later basis {T^} is built as the projective representa- 
tion of {Z/NZ®Z/NZ) in GL(iV,C) (see (jR4]l ). The corresponding structure constants (jBTTl) 
provides the Poisson-Lie brackets: 



(2.1) 



Let us introduce now the Lax matrix defined on the elliptic curve = C/(Z + rZ) with 
modular parameter r (Im(r) > 0): 



(2.2) 
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where = Z]^ (see (|B.3p ) and functions {fa{z — -^c)} form the basis in the space of sections 
r(Endy, St-) with simple poles at {zc}, c = 1, . . . ,n for the holomorphic vector bundle V of 
degree one associated with the principle GL(A^, C) -bundle over T,^. In fact the Lax matrix is 
fixed by the quasiperiodic properties with (IB.ip . ()B.2p : 

Liz + 1) = QL{z)Q-^, L{z + t) = kL{z)K-^ 

and residues YleSzaL{z) = S"*. 

The invariants of the Lax matrix generate commuting Hamiltonian^ 



1 " 

— = {H2M^ - ^c) - H^,cEi{z - z,)) - Ho, 

c=l 

where i?2,c = 2^ ((5''^)^) = | X] ^a^-a the Casimir functions corresponding to the or- 



bits Or and the Hamiltonians are: 



^1.'^ = Y.i^^MS")) = - E E Sl^SM^a - z,), (2.3) 

C byte 

= ^ E E S'~aS^pM - ^ E E S'-o.SMzt - Zc), (2.4) 

where we use the following notations: p{u}^) is defined in (jB.Sp . functions '^■y{z) and f-y{z) 
in (jB.lOp . (jB.lip . We also define the linear operators: 

p : Sa ^ SaP{iOa), 0ab '■ Sa — ^ Sa^a{Za - Zb), fab ■ Sa — > Safa{Za - Zb). 

In the following we also use Ei : S'a — s- SaEi{uJa)- Note that 

^*ab = -^ba (2.5) 

in the sense that {S^^abiS^)) = -{S^'PbaiS")) due to (|X6]) . Similarly, = ha, P* = P and 
El = -El. 

The commutativity of the Hamiltonians with respect to ()2.ip follows from the underlying 
linear r-matrix structure of the Belavin-Drinfel'd type: rf2'{z,w) = ^ ipa{z — w)Ta'S>T-a 

Note also that the Hamiltonians Hi^a are not independent: 

n 1 " 

e^i.^=45:e(^"^^<=(^^)) ^ 0- 

a=l a=l Cy^a 

The appropriate number of independent Hamiltonians is achieved by taking into account Hq 
and all higher Hamiltonians. 

Let us write down equations of motion with respect to the Hamiltonians ()2.3p . ()2.4p : 



dt^S- = {Hi,a,S'^} = -E['5^'^ac(5^)], (2.6) 

dt.S" = {Hi,a,S'} = [S',^ba{S'')], (2.7) 
dt.S'^ = {Ho,S'^} = [S^,p{S^)] - E[5^/ac(S^)]. (2.8) 



cj^a 



*Note that we use both the Eisenstein and the Weierstrass functions. They are simply related 
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Proposition 2.1. The equations of motion (j2.6p - (|2.8p can he presented in the Lax form ()1.6p 
with the Lax matrix L{z) defined in (|2.2p and M -matrices given as follows: 

Ma = S^^T^ifaiz - Za), (2.9) 

n 

Mo = -Y,Y. - ^b)- (2.10) 

Proof. The proof is direct. It is based on the usage of (|B.16p - ()B.20p . 

Let us prove identity ()B.20p which is the most nontrivial here. For a generic point w £ 'Sr 
consider m^{z,w) = (f^iz — w)(p^{w — Za)'- 

m'^ 1^1 ^^^^ _ ^^^(-^^(^ _ ^ Ei{w - Za) + Ei{ujy) - Ei{uj^ + z - Za)) 

= ip^{z - Za)iEi{z - w) + Ei{w - Za)) - f^{z - Za). 

Combining ()B.16p and (|B.19p which are implied to be known we have: 

ifl3{z - Zc)m'^ = {ipp{z - Zc)(p^{z - w))ip^{w - Za) 

ifisiw - Zc){^p+.y{z - w)ipj{w - Za)) + tp/B+jiz - Zc){iPj{Zc - w)ipj{w - Za)) 
llR16yR19|| ^^^^ _ Zc)(pp+y{z - Za)<f-p{w - Za) + tpp{w - Zc)ip^{z - Za)ipi3{z - w) 
+ ^P/3+-y{z - Zc)f^{Zc - Za){Ei{Zc - w) + Ei{w - Za) + Ei{u).f) - Ei{uj^ + Z^. - Za)) 

Lp^{z - Za)nfp - ip(S+^{z - Za)(pp{Za - Zc){Ei{w - Zc) + Ei{Za - w) + Ei{ujp) 

- Ei{ujj3 + Za- Zc)) + ^I3+^{Z - Zc)(p^{Zc - Za){Ei{Zc - w) + Ei{w - Za) + Ei{uj^) 

- Ei{uJ^ + Zc- Za)) 

= ip^{z - Za)rrff^ - ipi3+y{z - Zc)f^{Zc - Za) + ^/S+'yiz - Za) f jiiZa - Zc) 

+ {Ei{Zc - w) + Eii^W - Za)){iPl3+'y{z - Zc)^p^f{Zc - Za) + ifp+^fiz - Za)'.pj3{Za - Zc)) 

ip^{z - Za)m''p - ip(3+y{z - Zc)f^{Zc - Za) + ^13+yiz - Za)fl3{Za - Zc) 
+ {Ei{Zc -w)+Ei{w - Za))ip^{z - Za)ifj3{z - Zc). 

This ends the proof of ()B.20p . ■ 
2.2 Useful reformulation 

In this subsection we rewrite the equations of motion in a form which will be convenient for 1+1 
generalization. First, consider the following expressions for a = 1, . . . , n: 

r^(^-y(z - Za) ^ S'^if.yiZa - Zc) 
7Gr^ CT^a 

^7 S'^^'){Z - Zc){Ei{z - Za) + Ei{Za - Zc) + Ei{uj^) - Ei{z - Zc + UJ^)) 

= Ei{z - Za){L - Ma) + ^ McEiiZa " Zc) + Mq + T^S^f^{z - Za) 

c^a 7gr^ 

= Ei{z - Za)L + Y,McEi{za - Z,) + Mq - ^ T^S^^F^{z - Za). (2.11) 
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Then let us define new M-matrices in the fohowing way: 



Ma = ^ T^S!!^F^{z - Za)+ ^ T^7]l^(f^{z - Za), O = 1, . . . , 71, 

where 

ry"^ = T^S'^'p^iza - z^) = Y, M'iza) = Res,=,, (^—Liz)) . (2.12) 

From ()2.1ip we can see that the new M-matrices are the linear combinations of (j2.9p . ()2.10p : 



Ma = Ei{z- Za)L + McEiiZa " ^c) + Mq. 



Then the Lax equations yield 



a- L 



L^Y^cEiiza- Zc) + Mo 



J2 Ei{za - Zc)dt,L + dtoL 



and the equations of motion are: 

Y[S\E^{Za - Z,)^ac{S'^) - /ac(5^)] + [S^ p(S'^)] 

while for b ^ a: 

dij' = dt,S'Ei{za -Zb)+Y dt.S'Eiiza - Zc) + dt.S'' = [S\ pS""] 

c^a,b 

= [S', p{S') + Elizb - Za)^ba{Sn - fbaiS")] 

+ Y [S', {Elizb - Za) + EiiZa - ZcMbciS') - Ac(5')]. 
Cy^a,b 

Finally, we have 

di^S'^ = [S'',pS''] + Y[S'',Eac{S')], 

di^S' = [S',pS'] + [S\FUSn] + Y [S',M^ac{S''))] 

= YlS\M^aciS''))] + [S\FUSn] = [S\0baiv'n] + [S\FUS'')]. 

The corresponding Hamiltonians are obtained in the same way: 



b^c 
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= ^E^'5'f^^') + ^E^^"^-(^')) + ^ E iS^'MMS'))). (2.13) 

c CT^a b,c^a, bjtc 

The last one term equals: 

b,c^a, by^c b,cj^a cj^a 

= ^ E (^'^-(^a.(s^))) - ^ E<^'^(^')) + ^ E<^'^')^(^'^ - (2.14) 

b,c^a c^a c^a 

From (j2.13p , (|2.14p we conclude that the Hamiltonians for the reformulated version of the Gaudin 
model are of the form: 

Ha = ^{S^PS'^) + ^ Y.^S'^^aciS')) + ^Y. {S'^cai^abiS'))), a = 1, . . . , n 



c^a b,c^a 



or 



^- = W E^^'P^') + ^ E(^"^-(^')) + ^ E {S^^ca{^ab{S'))), a = 1, . . . , n. 

c c^a b,cj^a, 

b+c 

Two last forms of the Hamiltonians are differ by the constant 2^ ^ {S'^S'^)p{za — Zc). Let us 
summarize the obtained in results in 

Proposition 2.2. The dynamics of the Gaudin model produced by Hamiltonians 

Ha = ^{S^PS'') + ^ Y.(^^^ac{Sn) + ^ E {S'^PcaiMS'))) (2.15) 

cj^a b,cj^a 

is given by equations 

c^a 

and can he presented in the Lax form with L{z) from (j2.2p and 

Ma= Y Tf^i^-yi^ -Za)+ Y TyV^-yiz - Za), a = l,...,n. 

The Gaudin Hamiltonians ()2.3p and ()2.15p are simplified when written in terms of rj'"" ()2.12p : 

Ha = -{SS"), Ha = ^{S'^pS^) - ^iW)') + ^ Yi^^FaciS')). 

Cy^a 
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In the end of the section let us also give the rational "reformulated" version since it is more 
illuminating: 



Mc 1 ^ 1 ,^ V 



fa 



, Z Za Z Za Z Zq 



where 



Hamiltonians: 



2 



Ha 



'2N^\^Za-zJ N^^{Za-Z,y 

The later follows from simple evaluation: 



^ Za-Zc N ^fif Za- ZcZc- Zb N ^ {Za - Z^Y 



^ . (^'^ - ^«)(^^ - 2iV^^^ ^ ^ ^ AT 4> (z, - z,)2 ^ 2iV 4- (z, - z,)2 • 

The last one term is the analogue of the constant 2^ ^ {S^S^)p[za — Zc) in (j2.14p . The corre- 
sponding equations of motion are: 



{Za - ZbY Za- Zh^ Za- Zc (Za - Zb)"^ Za - Zb ' 



3 Field version 



3.1 1+1 s\{N,C) Gaudin model 

The general construction of the field version for the Hitchin systems was described in |19j . 
For our current purposes we only need to define the phase space. By analogy with mechanics 
let us consider a collection (direct product) of n orbits assigned to the marked points, i.e. let 
Res^=^^L(z) = S°'{x) be elements of the loop coalgebras si (A^, C) and x be a loop variable. We 
imply that the values of the invariants under the coadjoint action (or the eigenvalues of «S"^) are 
fixed. More over we assume for simplicity that the eigenvalues are C- numbers (independent of x). 
From the physical point of view it means that the magnetic momentum vector is normalized (as 
it is assumed in the Landau-Lifshitz model). The boundary conditions are chosen to be periodic. 
In summary, S"'{x) are si (A^, C)-valued periodic functions on a unit circle § : 5"(x+27r) = 5"(x) 
with eigenvalues {Xk,a, k = 1, . . . , N, a = 1, . . . ,n} fixed to be C-numbers: dxXa = 0. 
In the field case the Lax equations (jl.6p a replaced by the zero-curvature equations: 



dt^L-kdxMa = [L,Ma]. 



(3.1) 
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In fact, the numeration of Ma should include two type of indices as in (j2.3p : the first one type 
describes the number of the flow in the hierarchy and runs over 1, . . . , in 0+1 mechanics or 
1, . . . , oo in 1+1 field theory while the second one runs over 1, . . . , n in both cases and describes 
the assignment of the Hamiltonians to the marked points. In this paper we are not going to 
concern the whole hierarchy but only two first flows (as we did in 0+1 case). 

We will see that the first n flows of the hierarchy corresponds to the Gaudin Hamiltonians 
in the standard description ()2.3p supplemented by the momenta Pa along x while the second n 
flows naturally related to reformulated version (j2.15p 



Standard description Ha 



flows 



Reformulated version Ha 



ind 



flows 



Thus we do not use multi-index for times. It is sufficient to use ta and ia for our purposes 
and we keep these notations for the field version. 

It should be mentioned that the field generalization of the Lax pair into "L-A" pair satis- 
fying (j3.ip is nontrivial. The fact that the L-matrix (j2.2p is unchanged in the field version 
follows from the triviality of the moduli space of bundles of degree one. It is explained in [19j 
in detail. As a result we deal with the following Lax matrix: 

^ = E E T^S'^V^{z-z,). (3.2) 

The Ma-matrices for the first flow coincide with the mechanical versions either: 

Ma=Y, TjS;^j{z - Za). (3.3) 

Proposition 3.1. The zero- curvature equations ()3.ip with L from ()3.2p and Ma from (13. 3p are 
equivalent to the following equations: 

dt^s'^ - kd.s'^ = - J^[s^(^ae(s^)], 

dt^S'' = [S\ipUSn]- (3.4) 

The proof is the same as in the 0+1 case. As we will see below the Hamiltonian corresponding 
to Ma has the form 



na= (f dx{Pa + Ha{S{x))), 
J§1 



where ^gidxP^ is the shift operator in the loop algebra sl(A^, C): {^^i dx Pa{x), S''{y)} = 
SabdyS^{y) and Ha is defined as in (jl.Sp or (j2.3p . Thus the Hamiltonian describing equations 
(j3.4p has the form: 



-Ha = £ dx (^Pa-^Y^^ {S'^^aciS'))^ 



The phase space is a direct product of the symplectic orbits of the loop group SL(A^, C) with 
linear Poisson structure: 

{5^(x), 5^(y)} = 6abS{x - y)ca,fiS^^p{x), a,b = l,...,n. 

The second flows are of our main interest. 
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Proposition 3.2. The zero- curvature equations 

di^L-kd,Ma=[L,Ma] 

with L from (j3.2|) and 

Ma= Yl TjS^^ii^ -Za)+ J2 ^7<'/'7(^ " ^a), a = l,...,n, (3.5) 
where rf = r]'^ + Ar/", r]'"" = ^ Lpac{S^) are equivalent to the following equations: 

di^S'' = [S',^U^n] + [S\ha{S'')], 

-kd^S" = [S'' , Ar]"]. (3.6) 

The proof is also similar to the one given for the 0+1 case. Functions r]"" are not uniquely 
defined by equations —kdxS"" = [S'",A7/'^]. We fix this ambiguity by requiring r]"" — )• 77'" = 
^aciS") or A?7" in + 1 limit. As for the equation -kd^S'' = [5"^, Ary"] itself only 

some special cases were studied such as "vector" case [331 El] and "Grassmannian" case (special 
coadjoint orbits) [53] • For sl(2,C) case the answer is well known: Aif = —-^[S'^,S^]. 

3.2 2-site case and principal chiral model 

L-A pair for the principal chiral model was suggested in [29] (see also [231 [301 [361 EZ]). Consider 
the first flows of the Gaudin model (j3.4p with 2 sites or marked points (n = 2). It is convenient 
to start from the rational version: 

L = \ = Ml + M2. 

Z — Zl Z — Z2 

The corresponding M-matrix is known to be 

, , , , s s 

M = Ml- M2 = . 

Z — Zl Z — Z2 

Therefore the equations of motion are 

dtS^-kdxS^ = —[s\s% 

Zl - Z2 

dtS^ + kd,S^ = [S^ , . (3.7) 

Zl — Z2 

Then the Hamiltonian describing equations (j3.7p has a forir|^: 

{S^S'^) 



n = ni-n2= (p dxiPi-p^ 



Zl - Z2 

and the phase space is a direct product of two symplectic orbits of the loop group SL(A^, C) with 
the linear Poisson structure: 

{S^{x),S'p{y)} = 6ab6{x - y)c^,lsS^^+^ix), a, 6= 1,2. 



5 See Section 113] and g^Hl. 



12 



A.V. Zotov 



Remark 3.1. One can make a substitution = ^{Iq + li) and S"^ = ^(/q — ^i) to represent 
equations (|3.7p in its traditional form 

2 

dth - kdJo H [^1, ^o] = 0, 

Zl - Z2 

dth - kdJi = 

or change the coordinates {x,t) to "light-cone" coordinates ^ = ^i^-^, ry = *~^2 •" 
drtS' = —[S\S% 

Zl- Z2 

d,S^= ^ [S\S\ (3.8) 



Elliptic case. For L-A pair L = M1+M2 and M = M1-M2 with Ma = TaS^ipa{z-Za), 
a = 1,2 the equations ()3.4p yields {dt = dt^ — dt^)'- 

dtS^ -kd^s^ = -2[s\^i2{s^)i 

dtS^ + kd^S'' = 2[S\02i{S^)]. 
or by analogy with (jS.Sp 



dr^S^ = -2[S\^12{S^% 

d^S^ = 2[S\02i{S^)\. 

In sl(2,C) case this result was obtained by 1. Cherednik [30j. Here we see that the principal 
chiral model corresponds to the special (2-site) case of the first flows of 1 + 1 Gaudin model. 
It should be also mentioned that in [30] the equations for sl(2, C) case were obtained as a field 
version of XYZ model, i.e. from the second flow of 1-site Gaudin model (or sl(2, C) elliptic top). 
It may be explained as follows: consider stationary solutions S"" = S"'{ri) (or d^S"" = 0). Then 
fixing the ambiguity in solutions of the equation [S^ , ip2i{S^)] = as S*^ = —^^P2i{S^) we have 

dr,S' = [S\^l202liS^)] = [S\p{S% 

which is the equation of sl(A^, C) elliptic top (or 1-site elliptic Gaudin model) corresponding to 
the second flow H = {S^ p{S^)) . 

4 sl(2, C) 1+1 Gaudin models 

4.1 1+1 XXX Gaudin magnet: interacting Heisenberg models 

Let us consider the case Res2^-L(z) = S"* G sl(2,C) in detail. The linear Poisson~Lie structure 
in this case: 

3 

where 5^ are coefficients in the basis of Pauli matrices: 5"^ = ^ S^cja. 

a=l 

The Gaudin Hamiltonians are: 

_ _ ^ Oi Ji i- 02^2 "T '-'3 '-'3 

, Za Zc 
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while the Hamiltonians of the reformulated version are 
Since 

-^g _^ = _i ^ ('S"'>s''') _^ 

^Za-Z^ 2X1 2^ ^ Za- ZcZc- Zb 2X1 

~ 2 ^ {Za - Z,f ^ 2 (Ze-Za)(^c-26) ^ 

CT^a b,c^a\ bfc 

1^ 1 ^ (g^g^) , 



2^{Za-ZcY A iZa - Z^){Za - Zb) 2X1 

c^a b,cfa\ bfc 



then the corresponding equations of motion are: 

{Za - ZbY Za- Zh^^Za- Zc Za - Zfe (^a - Z^Y Za - Zb ' 

where 



,- = 5:^ + §5^ (4.4) 
Remark 4.1. (14. 4p differs from (j2.16p by -^5*" and the corresponding Hamiltonian (j4.3p differs 

a 

from (I2.17P by This difference does not follow from ansatz p.Sp but appears from the 

Hamiltonian description (see Section I4.3P . The corresponding Lax pair is given by L from (jl.2p 
and 



, Za Zc Z Za Z Za Z Za 

c=f:a 



1+1 version. Let S'^{x) G sl(2,C) be periodic sl(2, C)-valued functions on a circle S"^: 
S"'{x + 27r) = S°'{x) with eigenvalues {Aa} fixed to be C- numbers: dxXa = 0. The Poisson 
structure now is 



{S^{x),S'p{y)} = 2V^6'''ec,p^S^{x)6{x - y). (4.5) 
Consider 

Ma = 7 + 



(z - Za)^ Z- Za 

where 

= -IT2 [^"' Sl\ + 5] + = + 5^ = a.5^ (4.6) 
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Let us remark here that in (0+1) limit ry" = r(°' (j4.4p . Then the zero-curvature equation 

di^L-kd,Ma = [L,Ma] (4.7) 
reads as fohows 



E 



{Za-ZbY Za-Zb' 

These equation generahze the Heisenberg model which appears from 
case: 

dtS + ^[s,s^^] = o 

and described by the Hamiltonian 



n 



16A2 



dx {{d^sy 



(4.8) 
(4.9) 

in n = 1 (1-site) 



4.2 1+1 XYZ Gaudin magnet: interacting Landau— Lifshitz models 

By analogy with the previous section the Hamiltonians in 0+1 sl(2,C) case: 

= \{S^piS'^)) + ^ J](5"F(S^)) - i( (Y^^aciS'^)] ) + 
Consider now the following L-A pair: 

71 3 

H^) = YY1 S^f^o^ifaiz - Zc), 
c=l a=l 
3 

Ma{z) = ^ rj'^aaiPaiz - Za) + 5^CJq,97/3 (^ - Za)iP'y{z - Za), 



(4.10) 



where a, /3, 7 are different indices equivalent to 1, 2, 3 up to a cyclic permutation and (compare 
with gSDfl 



4A2 



(4.11) 



The zero curvature equation (I4.7P leads to equations of motion: 

5,-^5'^ - kd^v'^ = [S^ p(5«)] + MS')] - 0ca{[S', MS")]), 



drS' = [ipab{7ll,S'']+ipbama{S'),S-]) 



b\ aai 



(4.12) 



®The remark about the term ^S"' in the previous section is reasonable here as well. 
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and if (|4.1ip is a particular solution of the equation 



It is fixed if we require rj'^ ^ rj"- = ^ (pac{S'^) + -ttS"" in (0+1) limit. The proof follows 

from (ICl])-(Ia7l). In particular, from ([05]) it follows that ipcai[S^ , ^caiS"")]) = [Fca(5^), 5"]- 
Then (14.21) is written in the form close to (13.61): 



(4.13) 



dij'^ - kd.v" = [S^p(S'^)] + J][r?",^e,(5'=)] - [FUS''),S% 

di^S'' = [ipa,{7^%S'] + [S',FUS'')]. 



The last three terms in the first equation of (|3.6p vanish. It is due to ()C.3P that W A,B £ 
sl(2,C): 

Eii[A,B]) = [EiiA),B] + [A,EiiB)]. 
Note that in case n = 1 the first one equation in (|4.13p is the Landau-Lifshitz equation: 

dtS+^[S,S,,] = [S,p{S)] 
described by the Hamiltonian 

H = £ dx (^\{Sp{S)) + ^((9.5)2)) . 
4.3 Hamiltonian description 

The explicit form of the conserved quantities in terms of the fields are obtained by solving the 
Riccati equation. First, we make the gauge transformation (see for example [19^ I32j): 



kdx + 



L21 —Lii 



■01 

tp2 







kd:, + 



1 
T 



01 
1P2 







with 



T = L12L21 + Lj^ + kLu^^ - kd^Lu 



k'^dlLi2 ^ 3A;2 f d^Li2 



(4.14) 



L12 2 L12 4 y L12 

which leads to the Schrodinger equation: 

[-k'dl + T)^i = Q. 

Taking wave function in the form 0i = ■^==0 '^^ ^^^^ we come to the Riccati equation: 

kd^X + x" -T = 0. (4.15) 
The solution is obtained via local decompositions: 



Xa 



Z - Zn 



-Xa-1 + Xafi + {Z- Za)Xa,l + 



(4.16) 
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T=-. ^—.Ta,^2 + ^—Ta,-l+Ta,0 + ---, (4.17) 

[Z-Zay Z-Za 

L = -^—L"^-^ + L'^'O + {z- Za)L'''^ + ■■■ . (4.18) 

Z — Za 

Then (j4.15p gives: 

Xa -1 = -^a = Ta-2, 

Xafi = C^a -1 - kd^Xa-l), 

Xa,l = - ) , (4.19) 

2Xa,-l V 4ra _2 / 



As it was shown in [32] Xa.k (|4.19p are the densities of the conservation laws. We will use 
notation 

ha,k{x) = -XaXa,k-l 

for the densities and 



'Ha,k = f dx ha,k{x) 



for the Hamiltonians. The coefficients (14.18P of the decomposition of L-matrices in rational and 
elliptic cases: 



In what follows we sometimes omit for the simplicity the index a for L-matrix and its elements 
assmning decompositions (j4.16p - (|4.18p . Substituting (j4.18p into (j4.14p we get: 



r ,_r-lrO ,rO r-l.orO r-1, , -^n Ox^^u la t -1 



_ ^-1 r-l,r-lrl , 9 r 1 r-1 1 rO rO , rO rO 
-fa,0 — -'^12-'^21 "T ^12 ^21 -^-^ll^U ' ^12-^^21 -'^ll-'^ll 



+ -—1 ( LiidxL^2 + I^iidxLi2 — ^^irj — — kdxLii 

^12 V -t^l2 / 



9x^12 _|_ OxL^2 



2 L-|^2^ V -^12^ 



Let us summarize the obtained results. 
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Lemma 4.1. The density of the Hamiltonian Xa,o has the following form in terms of the de- 
composition of the L-matrix (|4.18|) : 

where Pa is the density of the Hamiltonian of the shift operator along x corresponding to the a*^ 
site {marked point) : 

[I dx Pa{x),S\y)\ = kSatdyS-, Pa{x) = - ^L^f' (x) ^f^^^ • (4.21) 

Proof. The nontrivial part of the proof is related to the shift operators Pa- Brackets (|4.5p have 
the following form in standard basis (here we omit index a as above): 

{L^iix),L^,Hy)} = 2L^i{x)6{x - y), {L^,\x) , L^,\y)} = -2L^l{x)5{x - y), 
{L-i(x), L-,\y)} = -4L-,\x)5ix - y). 

For example, let us verify that 4i dx | L^^ (x) , L2i{y)\ = —2dyL2i{y): 

L ^12 J 

f dx|Lr,Hx)MiM,L-/(y)| 

= I Ax {L-,l{x),L,^{y)f-^j^^ - i dxa.Lr/(x){lnLr2^(x),L2/(y)} 
= 2L2,Hy)^0^ + 4L^,Hy)^0^ = -25,L,"/(y). 

Li2 [y) ^12 iv) 

The later follows from the condition dyXl = = dy{{L^l {y))'^ + L^2iy)^2iiy))- The 
verification of (j4.2ip for other components can be performed in the same way. ■ 

Lemma 4.2. The density of the Hamiltonian Xa,i has the following form in terms of the de- 
composition of the L-matrix (|4.18p ; 

+ 2A:(L"'°a^.L'^'-iL'^'-i) - — ((9^L'^'-i)2). (4.22) 

Proof. Our purpose is to show that dx l.h.s. (|4.22p = <f^i dx r.h.s. (j4.22p . The proof is 
direct. It is based on the integration by parts. For example, 



dx 4X1 ( -y ^ + ^ (^)) - {^^^^^ - 



2 



- / dx -{L~'L~'' 
Jsi 2 



and 



2k (Bl^^ - d^L^A iL^2'4i + ^?2^2"i' + 2L?iL^ii) = 2k<[ dx {L°L-'L 
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Finally, we have the following densities of the Hamiltonians describing the first and the second 
flows: 

ha,l = -AaXa,0 = --^Ta-1 = Pa - - (L"'~^L°'°) , 

^a,2 = -KXa,l = {^a,0 - -^T^-^ = --(L°'°L"'°) 

These equalities are understood in a sense that dx l.h.s. (j4.23p = dx r.h.s. (|4.23|) . 

Theorem 4.1. The Hamiltonian densities (|4.23|) and the Poisson structure (j4.5p provides equa- 
tions dMD and (fiJ3|) . 

Proof. Substituting (j4.20p into (|4.23p we obtain explicit expressions for the Hamiltonians. In 
rational case: 

na,l = / dx ha,l{x) = I dx \Pa-\y] "^^^ | = / dx (P, + Ha), (4.24) 
T-ia,2 = f dx ha,2{x) 

= i'" |^^^-4A!5^:3^ + m!^^^^'^)^ 

and in the elliptic case: 

na,l = dx ha,l{x) = j^dx {Pa-\ J](5'^(^a,(5=)) j = dx (P, + F„), 
na,2 = dx ha,2{x) = j ^dx I \{S'^ piS'^)) + \ Y^{S'^F{S^)) " ^( I ^ '^-(^') ) 

where in Ha, Ha from ()4.ip . ()4.2p and (|4.10p S"^ = S"^(x) is assumed. The rest of the proof 
is simple. One should write equations of motion generated by the obtained Hamiltonians and 
the Poisson structure (|4.5p and verify that they coincide with equations (|4.8p and (|4.13p . This 
evaluation is direct. ■ 
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A Elliptic functions 

A.l Basic definitions and properties 

Most of definitions are borrowed from \38\ [39] , We assume that q = exp(27rir), where r is the 
modular parameter of the elliptic curve S,- which is realized as C/Tr, T,- = Z © Zr. 
The basic element is the theta function (e = exp 27rz) : 

t}{z\t) = ^(-l)"e ( -n(n + 1)t + nz\ 

ngZ ^ ^ 

oo 

= qh-'f{e'''' - e"^"") ^[{1 - g")(l - g"e2^"^)(l - r^'^e-^^"^). (A.l) 

n=l 

The Eisenstein functions 

Ei{z\t) = dzlog'&izlT), Ei{z\t) ~ - - 2r/iz, 

z 

where 

3 ^ 1 24 r?'(r) 



m=— 00 n.=— 00 ^ ' ' \ / 

and r/(r) = q^ Y\n>o(.^ ~ Dedekind function 

E2{z\t) = -d,Ei{z\T) = dllog^z\T), E2{z\t) ~ ^ + 2^. 

Z'^ 

The higher Eisenstein functions 

It is easy to see that the even-numbered functions are even and the odd-numbered ones are 
odd: 

E2k{-z) = E2k{z), E2k+i{-z) = -E2k+i{z). (A.2) 

Relation to the Weierstrass functions 

C{z,T)=Ei{z,T) + 2rii{T)z, (A.3) 
p{z,T) = E2{z,T)-2r]i{T). (A.4) 

The next important function is 

z) = u), (j){—u,—z) = —(j){u,z). (A. 6) 

It has a pole at 2: = and 

1 z 

4>{u, z) = - + Ei{u) + -(^f (-u) - p{u)) + ■■■ . 
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Let f{u,z) = du4>{u,z). Then 

f{u,z) = (l){u,z){Ei{u + z)-Ei{u)). 
Heat equation 

dr(t>{u,w) - -^dudu,(t){u,w) = 0. 



(A.7) 



Quasi-periodicity 






+ = -i?(z), 


i9{z + t) = 


-- -q~h~^''''^{z), 


Eiiz + l) = Ei{z), 


Ei{z + t) 


= Ei{z) - 2m, 


E2{z + l) = E2{z), 


E2{z + t) 


= E2{z), 


(f>{u, z + 1) = (f>{u, z), 


4>(.u,z + t] 




/(n, z + 1) = f{u,z), 


f{u,Z + T] 


= e-2™/(n,z)- 



(A.8) 



The Fay three-section formula: 

(j){ui,Zl)(p{u2, Z2) - (piui + U2, Zl)(j){u2, Z2 - Zl) - (j){ui + U2, Z2)(t){ui, Zl - Z2) = 0. (A. 9) 

Particular case of this formula is the Calogero functional equation 

(t){u, z)dy(j){v, z) - (j){v, z)du(l){u, z) = {E2{u) - E2{v))<p{u + v, z), (A.IO) 

(l){u, z)(l){-u, z) = E2{z) - E2{u), (A.ll) 

(p{z,ui)(f){z,U2) = (l){z,ui + U2){Ei{z) + Ei{ui) + Ei{u2) - Ei{z + ui + U2)). (A.12) 

Another important relation is 

Z — w)4){ui — V, z)(j){u2 + V,w) — (j){ui — U2 — V, Z — w)(j){u2 + V, z)(j){ui — V, w) 



= (j){ui,z)(j){u2,w)f{ui,U2,v), 

where 

f{ui,U2,v) = Ei{v) - -El (til - U2 -v) + Ei{ui -v)- Ei{u2 + v). 

Taking limit U2 ^ ^ in (lAT3]l we obtain: 

(j){v, z — w)(j){ui — V, z)(j){v, w) — (j){ui — v,z — w)(j){v, z)(j){ui — V, w) 
= (j){ui,z){E2{v) - E2{ui - v)), 

which is equivalent to (jA.lOp due to ()A.7p . 

Theta functions with characteristics. For a, 6 G Q by definition: 



(A.13) 



a 
b 



(z,t) = J^e((j + a)2- + (j + a)(z + 6)j . 



In particular, the function i) (jA.ip is a theta function with characteristics: 



i){x,t) 



1/2 
1/2 



(x,t). 
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Properties: 



a 
b 

a 
b 



{z + 1,t) = e{a)0 



a 
b 

./2 7" 



{z + a'r, r) = e ( -a' - - a'{z + b) 





a + a' 




b 



a + j 
b 



a 
b 



The following notations are used: 9 



a/2 
6/2 



(^) 



: Oabiz). Then Ouiz) = from (|AT]l . 



B Lie algebra sl(Ar, C) and elliptic functions 

Introduce the notation 
eN[z) = exp ( — ; 

and two matrices 

Q = diag(eAr(l), . . 
/ 1 



A 



Let 



1- 

0- 
V 1 • 

Fat = zj^^ = {Z/NZ © Z/NZ), 



e]\f{m 

\ 



1 

/ 



1 • • • 1 J 1 



(B.l) 



(B.2) 



N 



7(2) 



\ (0, 0) 



(B.3) 



be the two-dimensional lattices of orders A^^ and A^^ — 1 correspondingly. The generators of the 
lattice Ftv corresponding to the elliptic curve Hr (generated by the lattice Tt) are and ^. 

The matrices Q°'^A^^, a = (ai, 02) G Z]y generate a basis in the group GL(A^, C), while Q^^A'^^, 



a 



(ai,a2) e 



^(2) 
(2) 



generate a basis in the Lie algebra sl(A^, C). Consider the projective 



representation of Z]y in GL(A^, C) 



TaTb — ejv 



2 

a X b 



(a X 6 = 0162 — 0261)- 



The natural Killing form is 

{T^Tp) =^ Tr(rQT^) = iV(5o,«+^ mod r^v' 
It follows from (jB.5p that 



(B.4) 
(B.5) 

(B.6) 
(B.7) 
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where 



vr 

Sa,/3 = 2V-lsin — (a x /5) 



are the structure constants of sl(A^, C). They obey the following properties: 

Introduce the following constants on Z*^^-*: 
'71 + 72^' 



N 



Mu.,) = E, ( 21±^1 , E,{u:,) = ( 21±^^ , p(.,) = p ( 21±^ ) (B.8) 



N J ' " n N / ' °^ ° V N 
and the quasi-periodic functions on T,^ 

cp^(z) = ^i^^^±^,z), (B.9) 

ip^{z) = eAr(72z)(^^(z), (B.IO) 

/^(z) = eAr(72z)(9»0(n, 2:)|^_ 7i+72^ = (/?^(2;)(£'i(w^ + z) - £'i(a;^)), (B.ll) 

= - /^(z) = ip^{z){Ei{z) + Ei{uj^) - Ei{uj^ + z)). (B.12) 

Function f-yiz) is an element of the basis in the space of sections with a simple pole at z = of 
the bundle End(y) for the holomorphic vector bundle V of degree 1. It follows from (jA.Sp that 

ip^{z + l) = eNi^2)^-y{z), ip^{z + t) = eNi-^i)ip^{z), (B.13) 

f^{z + l) = eN{j2)fj{z), f^{z + T) = eN{-Ji)fj{z)-2Tn(fj{z), (B.14) 

F^{z + l) = eN{-f2)F^{z), F^{z + T) = eN{-li)F^{z). (B.15) 

Function F^{z) is the quasi-periodic and has the second order pole at 2; = 0. 
Let us write down the Fay-type formulae. It follows from ()A.9p - (|A.12p that 



Lp^{z - Za)^Pl3{z - Zc) ^= (pi3+.,{z - Za)fl3{Za - Zc) + iPl3+'y{z - Zc)(p..,{Zc- Za), (B.16) 

^p(.z)f^{z) - ip^{z)f(s{z) ^9 ip,3+'y{z){p{U(3) - piuJ^)), (B.17) 

ipa{z)ip-a(z) p{z) - p{uJa), (B.18) 
ipi3{z)ip^(z) ipi3+y{z){Ei{z) + Ei{uJi3) + Ei{uj^^) - Ei{z + UJ(S + UJ^f)), 

(/?^(zi)(^^(z2) ip^{zi + Z2){Ei{zi) + Ei{z2) + Ei{uj:y) - Ei{zi + Z2 + uj^)). (B.19) 

The last one identity can be also rewritten (using ()A.6P which is z) = — (^_^(z) and 
(H^) as follows: 

lPj{zi)lp-^{z2) = -ipjizi - Z2){Ei{zi) - Ei{z2) + Ei{uj.y) - Ei{zi - Z2 + a;^)). 
We also need the following relation: 

-iPjsiz - Zc)f'y{z - Za) + 937(2 - Za)fl3{z - Zc) 

= -ipi3+^{z - Zc)f^{Zc - Za) + ^p+'^iz - Za)fj3{Za - Zc). (B.20) 
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C Lie algebra sl(2, C) and elliptic functions 

For SL(2,C) instead of we use the basis of sigma-matrices 

(To = Id, cri = ro,i, 0-2 = ^1,1, o's = -Ti^o, 
{o-a} = {o-o,o-a}, a = 0,a, Q = 1,2,3, 

+ 2 ' 2 

The standard theta-functions with the characteristics are 

^0,0 = ^35 ^1,0 = ^2, ^0,1 = ^4, 01,1 = Gi- 
For Q = 1,2,3 and {ua} = {§, |} 



(/3q,(z) = e{zdrUJa)- 



Table 1. 



a 


(1,0) 


(0,1) 


(1,1) 




0-3 


0-1 


0"2 


half-periods 


Wi = i 


W2 = 5 


, , _ 1+r 
^3 - — 


(Pa{z) 


02(2)ei(O) 

62(0)01 (z) 


e4(2)9i(o) 

04(0)01(2) 


03(2)0i(O) 
03(0)01(2) 



In sl(2, C) case some more properties appear in addition to the previously listed. In what 
follows Q, /?, 7 are different indices equivalent to 1, 2, 3 up to a cyclic permutation. Then 



uja + ojp = uj^ mod r2 and ipa+piz) = (P'y{z), 

ip-aiz) = ^a{z), (fai-z) = -(faiz), (C.l) 

EiiiOa) = -27r^/^5,wa. (C.2) 
Indeed, from (IXl) and (IXSl) we have -Ei{^) = Ei{^ - t) = Ei{^) + 2^^^/^. Then 

Ei{uja) + Ei{uj^) = Ei{lj^ + oji3), (C.3) 

{ipa{z)f = p{z) - p{uJa). (C.4) 



For small z: 

1 z 

"Paiz) = -- -p{u;a) + ■■■ . 

From ()C.3P we also have 

F„(2;) = ipi3{z)ip^{z) = -d^ipaiz). (C.5) 
The Fay identity ()A.9p reads: 

ip.y{z - Za)ipi3{z - Zc) = ifaiz - Za)^Pl3{Za - Zc) - (Pa{z - Zc)^'y{Za - Zc). (C.6) 

Combining (jC.6p we may get: 

Lpp{z - Zc)(pi3{z - Za)(Paiz - Za) = ifisiz - Za)ip^{z - Za)^j3{Za - ^c) 

+ i^a{z - Zc)(Pa{Zc " Za)(p(3iZc - Za) - (Pa{z - Za)iPa{Zc - Za)^-yiZc - Za) (C.7) 

ip^{z - Zc)iPaiz - Za)^p^f{z - Za) = (f/siz - Za)^P'yiz - Za)ip-y{Za - Zc) 

+ ifaiz - Za)^a{Zc - Za)(Pl3{Za - Zc) - iPa{z - Zc)^aiZc - Za)^-yiZa - Zc). (C.8) 



24 



A.V. Zotov 



Acknowledgements 

Author is grateful to M.A. Olshanetsky for useful discussions and remarks. The work was 
supported by grants RFBR-09-02-00393, RFBR-09-01-92437-KEa, RFBR-09-01-93106-NCNILa, 
Russian President fund MK-1646. 2011.1 and to the Federal Agency for Science and Innovations 
of Russian Federation under contract 14.740.11.0347. 



References 

[1] Gaudin M., Diagonalisation dune classe d'Hamiltoniens de spin, | J. Physiqu^ 37 (1976), 1087-1098. 

Gaudin M., La fonction d'onde de Bethe, Masson, Paris, 1983 (in French); Mir, Moscow, 1987 (in Russian). 



[2] 
[3] 

K 



[5] 
[6] 

[7 
[9 

[lo; 



Baxter R.J., One-dimensional anisotropic Heisenberg chain, Pwn. Physics\ 70 (1972), 323-337. 

Sklyanin E.K., Takebe T ., Algebraic bethe ansatz for the XYZ Gaudin model, iPhys. Lett. A| 219 (1996), 
217-225, |q-alg/9601028| 



Zotov A., Elliptic linear problem for Calogero-Inozemtsev model and Painleve VI equation, 
\Lett. Math. Phys. 67 (2004), 153-165, |hep-th/0310260l 

Levin A., Olshanetsky M., Zotov A., Painleve VI, rigid tops and reflection equation, | Comm. Math. Phys.\ 
268 (2006), 67-103, |math.QA/0508058| " 

Levin A., Zotov A., On rational and elliptic forms of Painleve VI equation, Amer. Math. Soc. Transl. Ser. 2, 
Vol. 221, Amer. Math. Soc, Providence, RI, 173-184. 

Fuchs R., Uber lineare homogene Differentialgleichungen zweiterordnung mit im endlich gelegne wesentlich 
singularen Stellen, Math. Ann. 63 (1907), 301-323. 

Schlesinger L., Uber eine Klasse von Differentialsystemen beliebiger Ordnung mit festen kritischen Punkten, 
I J. Peine Angew. Math.\ 141 (1912), 96-145. 

Reyman A.G., Semenov-Tian-Shansky M.A., Lie algebras and Lax equations with spectral parameter on an 
elliptic curve. Zap. Nauchn. Semin. Leningr. Otd. Mat. Inst. Steklova 150 (1986), 104-118 (in Russian). 

Belavin A. A., Drinfel'd V.G., Solutions of the classical Yang-Baxter equation for simple Lie algebras, Funk- 
tsional. Anal, i Pnlozhen. 16 (1982), no. 3, 1-29 (English transl: \Funct. Anal. Appl.\ l& (1982), 159-180). 

Zotov A.V., Levin A.M.; Olshanetsky M.A., Chernyakov Yu.B., Quadratic algebras related to elliptic curves, 
\Theoret. and Math. Phys:\ 156 (2008), 1103-1122. larXiv: 0710. 1072] 

Sklyanin E.K., Separation of variables in the Gaudin model. Zap. Nauchn. Sem. Leningrad. Otdel. Mat. 
Inst. Steklov. (LOMI) 164 (1987), 151-169 (English transl.: J. Soviet Math. 47 (1989), 2473-2488). 
Sklyanin E.K., Takebe T., Separation of variables in the elliptic Gaudin model, [Comm. Math. Pfeys. | 204 
(1999), 17-38, solv-int/980700 8 

[11] Sklyanin E.K., Generating function of correlators in the SI2 Gaudin model, \Lett. M ath. Phys\ 47 (1999), 
275-292, solv-int/9708007', 

Takasaki K., Gaudin model, KZ equation and an isomonodromic problem on the tours, [Xett Math. Phys.\ 



44 (1998), 143-156, lhep-th/9711058 



Feigin B., Frenkel E., Reshetikhin N., Gaudin model, Bethe ansatz and critical level, [Comm. Math. Phys\ 
166 (1994), 27-62, hep-th/9402022 

Gould M.D., Zhang Y.-Z., Zhao S.-Y., Elliptic Gaudin models and elliptic KZ equations, \Nuclear Phys. -B| 
630 (2002), 492-508, InhmST/OTToMsl, 



Chernyakov Yu., Levin A., Olshanetsky M., Zotov A., Elliptic Schlesinger system and Painleve VI, 
,7. Phy sTA: Math. Gen.| 39 (2006), 12083-12101, | nlin.SI/0602043| 

[12] Khesin B., Levin A., Olshanetsky M., Bihamiltonian structures and quadratic algebras in hydrodynamics 
and on non-commutative torus, \Corrim . Math. Phys. 250 (2004), 581-612, nhn.SI/0309017 

[13] Petrera M., Suris Yu.B., An integrable discretization of the rational su(2) Gaudin model and related sys- 
tems, I Comm . ' Math. Phys. 283 (2008), 227-253, arXiv:O707.4O88 

Ragnisco O., ZuUo F., Backlund transformations for the trigonometric Gaudin magnet, SIGMA 6 (2010), 
012, 6 pages, larXiv:0912.2456l 

Veselov A. P., What is an integrable mapping?, in What is Integrability?, Editor V.E. Zakharov, Springer 
Ser. Nonlinear Dynam., Springer, Berlin, 1991, 251-272. 

Hone A.N.W., Kuznetsov V.B., Ragnisco O., Backlund transformations for the sl(2) Gaudin magnet, 
J. Phys. A: Math. Gen. 34 (2001), 2477-2490, |nlin.SI/0007041, 



1+1 Gaudin Model 



25 



[14] Talalaev D.V., The quantum Gaudin system, Funktsional. Anal, i Prilozhen. 40 (2006), no. 1, 86-91 (Eng- 
lish transL: \Funct. Anal. Appl. 40 (2006), 73-77). 

Rubtsov v., Silantyev A., Talalaev D., Manin matrices, quantum elliptic commutative families and charac- 
teristic polynomial of eUiptic Gaudin model. C^T^MJI 5 (2009), 110, 22 pages. larXiv:0908.4064l 

[15] Frenkel E., Affine algebras, Langlands duality and Bethe ansatz, |q-alg/9506003 

Chervov A., Talalaev D., Quantum spectral curves, quantum integrable systems and the geometric Langlands 
correspondence, |hep-th/ 0604128 

Teschner J., Quantization of the Hitchin moduli spaces, Liouville theory, and the geometric Langlands 
correspondence I, arXiv: 1005.2846 

[16] Enriquez B., Rubtsov V., Hitchin systems, higher Gaudin operators and _R-matrices, Math. Res. Lett. 3 
(1996), 343-357, [aig-geom/95030l0| 

[17] Hitchin N., Stable bundles and integrable svstems. rPwfce Math.~I\ 54: (1987), 91-114. 

[18] Nekrasov N., Holomorphic bundles and many-body systems, \Comm. Math. Phys.\ 180 (1996), 587-603, 
|hep-th/9503157| 

Gorsky A., Nekrasov N., Elliptic Calogero-Moser system from two dimensional current algebra, 
,hep-th/9401021 , 

[19] Levin A., Olshanetsky M., Zotov A., Hitchin systems - symplectic Hecke correspondence and two- 
dimensional version, II Comm. Math. Phys.\ 236 (2003), 93-133, [nIIn.SI/0110045| 

[20] Petrera M., Ragnisco O., From su(2) gaudin models to integrable tops. \SIGMA\ 3 (2007), 058, 14 pages, 
|math-ph/0703044 . 

Levin A., Zotov A., An integrable system of interacting elliptic tops, Teoret. Mat. Fiz. 146 (2006), 55-64 
(English transL: | Theoret. and Math. Phya] 146 (2006), 45-52). 

[21] Kuznetsov V.B., Isomorphism of the n-dimensional Neumann system and the n-site Gaudin magnet, Funk- 
tsional. Anal, i Prilozhen. 26 (1992), no. 4, 88-90 (English transL: \Funct. Anal. Appl.\ 26 (1992), 302-304). 

[22] Zakharov V.E., Shabat A.B., A scheme for integrating the nonlinear equations of mathematical physics 
by the method of the inverse scattering problem. I, Funktsional. Anal, i Prilozhen. 8 (1974), no. 3, 43-53 
(English transL: \F \inct. Anal Appl. 8 (1974), 226-235). 

[23] Faddeev L., Takhtajan L., Hamiltonian approach to solitons theory, Nauka, Moscow, 1986 (in Russian). 

[24] Dubrovin B.A., Krichever LM., Novikov S.P., Integrable systems. I, Current Problems in Mathematics. 
Fundamental Directions, Vol. 4, Itogi Nauki i Tekhniki, Akad. Nauk SSSR, Vsesoyuz. Inst. Nauchn. i Tekhn. 
Inform., Moscow, 1985, 179-284 (in Russian). 

[25] Sokolov V.V., Schabat A.B., Classification of integrable evolution equations, Soviet Sci. Rev. Sect. C Math. 
Phys. Rev., Vol. 4, Harwood Academic Publ., Chur, 1984, 221-280. 

Mikhailov A.V., Schabat A.B., Yamilov R.I., The symmetry approach to classification of nonlinear equa- 
tions. Complete list of integrable systems, Uspekhi Mat. Nauk 42 (1987), no. 4, 3-53 (English transL: 
\Russian Math. Survey 8^ 42 (1987), no. 4, 1-63). 

Fokas A.S., Symmetries and integrability. Stud. Appl. Math. 77 (1987), 253-299. 

[26] Landau L., Lifshitz E., On the theory of the dispersion of magnetic permeability in ferromagnetic bodies, 
Phys. Zeitsch. der Sow. 8 (1935), 153-169. 

[27] Sklyanin E., On complete integrability of the Landau-Lifshitz equation, Preprint LOMI E-3-79, 1979. 

Borovik A.E., Robuk V.N., Linear pseudopotentials and conservation laws for the Landau-Lifshits equation 
describing the nonlinear dynamics of a ferromagnet with uniaxial anisotropy, Teoret. Mat. Fiz. 46 (1981), 
371-381 (English transL: Theoret. and Math. Phy s. 46 (1981), 242-248). 

[28] Baxter R.J., Eight- vertex model in lattice statistics and one-dimensional anisotropic Heisenberg chain. 
I. Some fundamental eigenvectors, Ann. Physics 76 (1973), 1-24. 

Takhtajan L.A., Faddeev L.D., The quantum method of the inverse problem and the Heisenberg XYZ model, 
Uspekhi Mat. Nauk 34 (1979), no. 5, 13-63 (English transL: j Russian Math. Surveys\ 34 (1979), no. 5, 11-68). 

[29] Zakharov V.E., Mikhailov A.V., Relativistically invariant two-dimensional models of field theory which are 
integrable by means of the inverse scattering problem method, Soviet Phys. JETP 74 (1978), 1953-1973. 
Pohlmeyer K., Integrable Hamiltonian systems and interactions through quadratic constraints. 
Com m. Math. Phys.\ 46 (1976), 207-221. 

[30] Cherednik I.V., Local conservation laws of principal chiral fields (d = 1), Teoret. Mat. Fiz. 38 (1979), 179- 
185 (English transL: \ Theoret. and Math. Phys] 38 (1979), 120-124). 

Cherednik I.V., Relativistically invariant quasiclassical limits of integrable two-dimensional quantum models, 
Teoret. Mat. Fiz. 47 (1981), 225-229 (English transL: ^Theoret. and Math. Phys. ,47 (1981), 422-425). 



26 



A.V. Zotov 



[31] Krichever I., Vector bundles and Lax equations on algebraic curves, I Comm. Math. Phy s]229 (2002), 229- 
269, |hep-th/0108110[ ^ 

[32] Dubrovin B.A., Matveev V.B., Novikov S.P., Non-linear equations of Korteweg-de Vries type, finite-zone 
linear operators, and Abelian varieties, Uspekhi Mat. Nauk 31 (1976), no. 1, 55-136 (English transL: 
\Russian Math. Surveys\ 31 (1976), no. 1, 59-146). 

[33] Golubchik I.Z., Sokolov V.V., Multicomponent generalization of the hierarchy of the Landau-Lifshitz equa- 
tion, Teoret. Mat. Fiz. 124 (2000), 62-71 (English transL: Theoret. and Math. Phys. 124 (2000), 909-917). 

[34] Akhmetshin A. A., Krichever I.M., Volvovski Yu.S., Elliptic families of solutions of the Kadomtsev- 
Petviashvili equation, and the field analogue of the elliptic Calogero-Moser system, Funktsional. Anal, i 
Prilozhen. 36 (2002), no. 4, 1-17 (English transL; Funct. Anal. Appl.^ ^Q (2002), 253-266), ,hep-th/0203192i 

[35] Skrypnik T., 'Doubled' generalized Landau-Lifshiz hierarchies and special quasigraded Lie algebras, 
I J. Phys. A: Math. Gen. 1 37 (2004), 7755-7768. 

Skrypnik T., Quasigraded Lie algebras and matrix generalization of Landau-Lifshitz equation, in Pro- 
ceedinds of Fifth International Conference "Symmetry in Nonlinear Mathematical Physics" (June 23-29, 
2003, Kyiv), Editors A.G. Nikitin, V.M. Boyko, R.O. Popovych and LA. Yehorchenko, Proceedings of Insti- 
tute of Mathematics, Kyiv 50 (2004), Part 1, 462-469. 

[36] Holod P.I., The hidden symmetry of the Landau-Lifshits equation, the hierarchy of higher equations 
and a dual equation for an asymmetric chiral field, Teoret. Mat. Fiz. 70 (1987), 18-29 (English transl.: 
I Theoret. and Math. Phys . 70 (1987), 11-19). 

[37] Orlov A.Yu., A'^-soliton solution of chiral fields on Grassmann manifolds (a model), Teoret. Mat. Fiz. 61 
(1984), 214-225 (English transl.: ITfeeoref. and Math. Phys.\ Ql (1984), 1099-1107). 

[38] Weyl A., Elliptic functions according to Eisenstein and Kronecker, Ergebmsse der Mathematik und ihrer 
Grenzgebiete, Band 88, Springer- Verlag, Berlin - New York, 1976. 

[39] Mumford D., Tata lectures on theta. I, Progress in Mathematics, Vol. 28, Birkhauser Boston, Inc., Boston, 
MA, 1983. 

Mumford D., Tata lectures on theta. II. Jacobian theta functions and differential equations. Progress m 
Mathematics, Vol. 43, Birkhauser Boston, Inc., Boston, MA, 1984. 



